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Abstract: We begin by studying a very simple Hamiltonian for Maxwell’s equations that 
has no gauge fields and is made entirely of the electromagnetic fields. We then show that 
this theory cannot be quantized. We also show that no other such simple theory that only 
involves the electromagnetic fields can be quantized. This gives further evidence for the 
important role of gauge fields in QED. 
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1 Introduction 


The standard variational approach to Maxwell’s equations is based on the action 

S[A] = I d^x, (1.1) 

where the metric signature is-h ++, and = d^A^ — d^A^ is the held strength 

tensor. There are various physically equivalent approaches to quantization. One of the 
most convenient ones is that of constrained systems [1]. It begins by dehning the density 
of momentum conjugate to A^: 

:= or TTi = -Ei. (1.2) 

The Poisson brackets of Ai and vr^ obey the standard formula for canonical variables, from 
which the Poisson brackets of Ei and Bj are to be computed: 

{^i(x), Aj(x')} = {7ri(x),7rj(x')} = 0, {Aj(x), 7rj(x')} = 5ijS^{x- x), 

{Ei{x),Ej{x)} = {Bi{x),Bj{x)} = 0 , {Ei{:>z),Bj{x)} = {-7ri(x),ejkidkAi{x)} 

= eijfc^<^^(x-x'). 

These Poisson brackets can then be promoted to commutators to quantize the theory. 

If we forget about the sources for a moment and set J = 0, we see that the action is 
a function of the helds E alone. Nonetheless, varying S with respect to F does not yield 
the correct equations of motion. It seems that A has a fundamental role without which 
Hamiltonian, action, and quantization cannot be defined. Of course, gauge symmetry has 
proved to be a cornerstone of particle physics and there is no doubt that it has wonderful 
theoretical and experimental consequences. But the question here is whether the quantiza¬ 
tion or at least a Hamiltonian formulation of Maxwell’s equations without resort to gauge 
fields is possible. Since ultimately, the electric and magnetic fields—rather than gauge 
fields—are the physical variables, this is a well-motivated question. 

There is a conventional argument for why quantization necessitates gauge invariance 
[2]. It is based on the particular nature of massless spin-1 particles, and leads to the 
conclusion that the causal field operator made out of their creation/annihilation operators 
has to obey gauge symmetry. Although this argument is compelling in many respects, it 
doesn’t forbid construction of the Hamiltonian out of the electromagnetic field F. The aim 
of this paper is to address this question. But we will not take the particle theoretic point 
of view where the mass and spin of the particle are crucial. Indeed, we solely concentrate 
on the Hamiltonian formulation of Maxwell’s equations and its quantization (which we call 
QED). 

The rest of this paper is organized as follows. In Section 2 we derive a Hamiltonian 
for the classical electromagnetism in which gauge fields play no role. In Section 3 we 
explain the trouble with the quantization of this model and its modifications. Finally, we 
summarize and conclude in Section 4. 
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2 The Gauge-less Formulation 

There are two of Maxwell’s equations, 


dE ^ ^ ^ 


V • E = p, 
V-B = 0, 


( 2 . 1 ) 


that contain time-derivatives and are reminiscent of Hamilton’s equations, as if E and B 
are canonical coordinates. They suggest that we may find a Hamiltonian, exclusively in 
terms of E and B, whose equations of motion are Maxwell’s equations. 

This motivates us to try to construct a gauge-less Hamiltonian of the form Hgi = 
f 9{gi(E, B)d^x and treat EjS and HjS as the generalized coordinates and momentum den¬ 
sities, i.e., 

{Ei(x),Bj(x')}gi = £~^6ijd^(x - x'), (2.2) 

where £ is an irrelevant length scale required on the basis of dimensional analysis. We note 
in passing that this is in sharp contrast with Eqs. (1.3); we will get back to this point later. 
We further demand that Tigi involve at most one spatial derivative and be quadratic in 
its arguments, so that the equations of motion are linear in EjS and HjS with first order 
spatial derivatives. With these restrictions, a generic Bgi reads 


£ ^T~Lg\ — EittijEj BibijBj EiCijBj -|- diEi -|- CiBi, (2.3) 

where Uij = a® -|- a^^ldk and (similarly for b and c), while di = df'^ + d^^^dj 

(similarly for e). 

To get the desired set of equations of motion in the left column of (2.1), we need 

= {-Fij-f^gljgh = {Bi,-f^gljgh (2.4) 

which translates to: 


eijkd.B, -.h = h,B, + B,bf - 4 {B,b%) + Ejcf - 4 + ef - 4eg\ 


This implies that 


jEj 


(e, 

41 ) - 

CijBj 

-((<“> +a, 








(2.5) 

Cij 

= di = 0, 


g(0) _ n (1) _ 




II 


''S ■ 

- = 

^ikj ? 

(2.6) 


= 9ka% 


_^(1) , „(1) _ 
^ijk ^ '^jik — 

^ikj • 



; + hgi with 

a'" 

%k 

- ^jik 

(similarly for 6); and the Hamiltonian 


density becomes: 


1 


r%i = ^E.VxE + ^B.VxB-J.B + ^4 [d^lE.E, + b\;>B,B,]+d, {e\;>B,) . (2.7) 




ijk 


ijk 


,( 1 ); 
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Of course, we should ignore the total derivatives—as we did so in the first place when 
writing Eq. (2.3)— since they disappear in our final result for the Hamiltoniand 

QE-VxE+^B-VxB-J-B^d^x. (2.8) 

Note that p and J are sources that have to satisfy the continuity equation, dp/dt+V-J = 0. 
(For the purposes of this paper, the sources can be considered non-dynamical. But one 
can associate their dynamics with electrically charged particles if those degrees of freedom 
are taken into account.) 

So far our Hamiltonian produces six equations out of the eight Maxwell’s equations. 
We must supplement the remaining two Maxwell’s equations (the right column of (2.1)) 
as constraints on initial conditions and ensure that they remain satisfied at all subsequent 
times. This requires: 

dp 

{diEi, Hgilgi = —, {diBi, Hgilgi = 0. (2.9) 

But these are trivially satished as a consequence of the original six equations (the left 
column of (2.1)) and as a result of the continuity equation. So once the initial conditions 
are chosen to satisfy the constraints, they are guaranteed to remain satisfied under the 
Hamiltonian evolution. This shouldn’t undermine the significance of the constraints; if it 
were not for them, the charge density would not appear anywhere in the solution, since H^i 
is independent of p. Finally, we note that the two constraints V • E — /? = 0 and V • B = 0 
are second-class, so it is not possible to obtain a more general motion by adding them to 

The reader may wonder how our Hamiltonian Hg\ is related to the standard Hamilto¬ 
nian,^ 

Hstd = j - J • d^x, (2.10) 

and energy of the electromagnetic held, 

£ = ^j{E‘^ + B‘^)d^x. (2.11) 

Of course, the gauge-less formulation has to yield the standard result that when there is no 
dissipation, S must be a conserved quantity. Indeed, using the gauge-less Poisson brackets 
in Eq. (2.2), it is straightforward to show that 

7 ^ p 

— = {£, Hgilgi = - y [V • (E X B) + J • E] d^x. (2.12) 

similar expression appears, e.g., in Ref. [3] and is based on the Lagrangian density C = — jF* • 
{idF/dt -I- V X F) -I- c.c. introduced in Ref. [4], in which F = B + iE. However, these works begin by 
assuming E and B as independent coordinates in a 6 -I- 6-dimensional phase space, which is different from 
our 3 -I- 3-dimensional space. Their action differs from the one we derive below in Eq. (2.16), too. 

^The most general Hamiltonian in the standard theory is more complicated. Here Hstd governs the time 
evolution in the temporal gauge (A° = 0). 
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Thus we have reproduced the integrated form of Poynting’s theorem d£/dt + / V • Sd^x = 
— f J ■ Eid^x. In a similar manner, using the standard Poisson brackets in Eq. (1.3), one 
obtains: 

^ ^ + {^gi, ^stdlstd = -^ / j • (2-13) 

This confirms that {Hgi, ffstdjstd = 0 as it should, since the time evolution of Hgi can be 
obtained by either i/gtd or by Hgi itself.^ 

Since we are going to use canonical quantization, we will not need a Lagrangian. But 
to complete this section, we make a brief comment. It is straightforward to obtain the 
Lagrangian density: 

£gi = E -B • V X B - • V X E. (2.14) 


Of course, £gi has to be considered as a function of E and E, which can be done, using the 
Helmholtz theorem, by writing B in terms of E: 


B(x) = V X 


E(xO + J(xO 3 
47r|x — x'l 


(E +J),,/ X (x-xQ 3 
47r|x — x'P 


(2.15) 


We also need to impose the constraint V • E = p by a Lagrange multiplier u(x) in the 
Lagrangian (there is no need to impose V-B = 0, since it is already satisfied by Eq. (2.15)— 
in the Lagrangian formalism B is defined by Eq. (2.15)). So finally the action corresponding 
to Hgi reads: 


5gi[E] = 


1 

Stt 

4 


X — X 
|X-X'|3 


(E + J)x X (E + J)^ 


d^x'd^x 


(E • V X E)^d^x + / m(x)(V • E - p)xd^x. 


(2.16) 


It is easy to see that 5Sg\ = 0 implies V^n = 0. If the fields vanish at infinity, 
the constraint is satisfied and we can have n = 0 there. = 0 then implies that tt = 0 
elsewhere. Thus we can set n = 0 in all equations of motion (after performing the variation) 
to obtain d^/dt = —V x E. Gauss’s law is imposed as a constraint, and the remaining 
two of Maxwell’s equations follow from the definition of the conjugate momentum in the 
action. Note how this differs from the standard formulation in which the sourced equations 
follow from the action, while the sourceless ones are automatic. Also note that we have 
suppressed the time dependence in Eqs. (2.15) and (2.16): all quantities are evaluated at 
the same time x^ = x^^. Therefore, we observe that although the action is local in time, 
it is very non-local in space. In contrast, the Hamiltonian Hg\ is local in both time and 
space. 

There is another aspect of the gauge-less formulation which is more important than the 
non-locality of action, and that is the lack of Lorentz invariance. In particular, unlike the 

®Note that H^\ cannot be used to generate the time evolution of -ffatd, since the latter involves gauge 
fields, whose Poisson brackets are undefined in the gauge-less formulation. 

^By “local in space” we mean the standard terminology that the second and higher spatial derivatives 
of the fields do not appear in T-i. 
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standard Poisson brackets in Eq. (1.3), our Poisson brackets in Eq. (2.2) are not Lorentz 
invariant. Despite this, the gauge-less formulation leads to Maxwell’s equations which are 
clearly Lorentz invariant. As long as we are only concerned about the equations of motion, 
none of these problems is an issue. However, as we see in the next section, Lorentz violation 
does obstruct the quantization process of the gauge-less formulation. 

3 Quantization 

As we remarked before, the fact that the Poisson brackets in Eq. (2.2) are not covariant is 
not an issue in the classical theory, since all that matters there is the equation of motion. 
When we quantize the theory, however, the commutator of two observables has physical 
meaning; it is related to the uncertainty after all. Let us inspect in more details this 
distinction between the classical and quantum theory. 

Quite generally, let and pa be the generalized coordinates and momenta describing 
a given system in one reference frame, and let (^°' and describe the same system in 
a second frame. Of course, we must have symplectic structures {,} and {,}' on both 
sets of coordinates satisfying {q°‘,Pb} = and {q'°‘-,p'i^' = 5^. Eurthermore, the primed 
coordinates are related to the unprimed ones by Lorentz transformations and we write 
q\q,p) and p'{q,p) to express this relationship. In most applications, the primed and 
unprimed Poisson brackets are the same, so that = {q''^{q,p),p'i,{q,p)}- But 

this is not a necessity in a classical theory. Our formulation of Maxwell’s equations in 
the previous section is one such example, where / {E'(£', H), (E, H)}. To be 

specific, under a boost in the x-direction, we have = 7 (E 2 —uEs) and Eg = ^{E 3 +vB 2 ), 
which implies {E 2 (x), E 3 (y)} = 2u7^(5^(x—y) ^ 0. Nevertheless, the gauge-less description 
can successfully reproduce Maxwell’s equations. 

The commutator of two observables is not as flexible as their Poisson bracket. The 
reason lies in the fact that probability is invariant and so the Hilbert spaces of the two 
observers must be unitarily related. In practice, it is common to use a single Hilbert space 
for both observers and employ identical operators to describe physical observables. Then 
we relate the state vector as seen by one observer by a unitary transformation to the 
state vector as seen by the other observer. Thus, there is only one kind of commutator; 
symbolically: [, ] = [, ]'. In particular, for the canonical operators and p^, we must have 
W'^{q,p),Phiq,p)] = [q'',Pb]- 

We are thus led to conclude that the gauge-less formulation of the previous section 
cannot be quantized since its Poisson brackets are not covariant under Lorentz transforma¬ 
tions. The natural question is whether it is possible to write down a set of covariant Poisson 
brackets involving the physical fields (rather than gauge fields). We now try to answer this 
question under the assumption that the canonical coordinates of interest are linear combi¬ 
nations of the electric and magnetic fields, to allow for the pairing of coordinates and their 
conjugate momenta to be different from that of Eq. (2.2). 

Specifically, let us assume that there are three pairs of canonically conjugate fields, 
collectively denoted by q{x)s and p^{x)s, which are linear combinations of the six nonzero 
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components of The qs and p^s must be organized as elements of some tensor in 

order for the commutators to be covariant. Therefore, the qs and p^s are linearly related 
to T^jyS with coefficients that are invertible tensors. We take the coefficients appearing 
in these linear combinations to be independent of x, to respect translation invariance. 
Inverting these relations, we can write F^uS as linear combinations of qs and ps. Since the 
commutators of the latter are given by Eq. (A. 7), the commutator of the former must be 
of the form; 

[Ff,^{x),Fpa{0)] = iAp^pa{x)A{x) + iBl^p„{x)drA{x). (3.1) 

The operator B cannot be identically zero, since [g,p]s are linear in [F, E]s after all. In 
fact, by virtue of the constraint on the coefficient of dpA in Eq. (A. 7), B must remain a 
nonzero tensor as x approaches the origin. Eurthermore, this nonzero tensor cannot depend 
on the direction of x (for example, it cannot be because we took the coefficients 

relating qs and ppS to F^j^s to be x-independent. Altogether, these imply that the tensorial 
structure of BJ^^^p^ must be made up entirely of the invariant tensors, namely, rjpy and 
This brings about a contradiction, since B has an odd number of indices and cannot 
be possibly made up of the even-ranked tensors p and e.® 

We conclude that the electric and magnetic helds themselves (or linear combinations of 
them) cannot provide a Lorentz invariant set of canonical variables. Note how the standard 
formulation avoids this conclusion: Since the electromagnetic fields are derivatives of the 
canonical variables (i.e., the gauge fields), their commutator contains two derivatives of A 
(equivalently, one derivative of <5^) and B is identically zero. There is no obstruction in 
constructing the tensor A^^p^j in Eq. (3.1) out of the metric and x^ (equivalently dpA). 
Indeed, for the free Maxwell theory, we have [Fpu{x), Fpc^{x')] = z9[^7yj,][o-9p]'A(x — x'). 

4 Conclusions 

We saw that a Hamiltonian formulation of Maxwell’s equations without resort to gauge 
fields is possible. This gauge-less classical theory has a local Hamiltonian, but its action 
is non-local in space. Although, the theory ultimately yields Lorentz invariant equations 
of motion, the action and the canonical structure are not preserved under Lorentz trans¬ 
formations. This is not a problem in a classical theory, but it does affect the quantum 
theory. 

We then showed that the gauge-less model cannot be modihed to have covariant com¬ 
mutators. We did so by considering all possible ways of writing three canonical variables 
and their conjugate momenta as linear combinations of the electric and magnetic fields. It 
turns out that this is inconsistent with what we expect from the commutators of electro¬ 
magnetic fields. Of course, more complicated relationships between the canonical variables 

^Useful properties of commutators of quantum fields in the notation appropriate for this discussion are 
summarized in Appendix A. 

®We have also assumed that there is no preferred 4-vector in the theory. We can turn off the sources 
for the sake of this discussion, since if—as we showed—the theory doesn’t work in the absence of the 
sources, it won’t work in their presence either. 
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and the electromagnetic fields can be studied. But we restricted our attention to the sim¬ 
plest possibility where the electromagnetic fields themselves are the building blocks of the 
canonical formalism. 
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A Commutators of Quantum Fields 

In this appendix we review some basic properties of commutators of generic (free and non- 
free) relativistic quantum fields. Since we have undertaken a Hamiltonian approach, we 
use the notation q and p for canonical variables even when they are field operators. 

Let n = (1, 0) be the 4-velocity of an observer, and x be the proper spatial coordinates 
on the surface t = 0. According to this observer, the canonical fields obey 

[gn(x),gn(x')] = [pn(x),Pn(x')] =0, (x),p„(x')] = (x - x')I. (A.l) 

For a scalar field, the Heisenberg field operator q{x) is obtained by evolving ^n(x) with the 
Hamiltonian Hn{qn,Pn) to any future or past time. It encapsulates all of the information 
in ( 7 n(x)s for all timelike n and all spacelike x. In particular, [g(x),g(0)] = 0 for spacelike 


X. For timelike x, this provides even more information and we have 

[g(x),g(0)]=zC('?'?)(x)A(x), (A.2) 

where the operator contains this extra information (for a free field of mass m, it is 

equal to the identity operator). 

A(x) = , [ sgn{k^)S{k‘^ + d'^k (A.3) 

(27r)^ J 

is Schwinger’s A function with the following properties: (i) A(—x) = — A(x), so that 
A(0) = 0; (ii) A(x) vanishes for spacelike x; and (iii) for any spacelike surface S passing 
through the origin and having normal vector 

[ d^/^{x)da^ = -1. (A.4) 

Jt. 

Unlike the g„(x)s which give rise to the n-independent q{x), the Pn{^) do not combine 
into a scalar operator p{x), rather into a 4-vector Pfi{x) such that 

[q{x),p^m = -iC^^P\x)d^^{x), C(''P)(0) = i. (A.5) 

Using Eq. (A.4), the desired commutation relations in (A.l) are obtained for pn = n^p^. 
Again the operator carries additional information (for a free field, p^ = d^q and 

C(<iP) = i). 
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The commutator of the momenta at spacelike separation must vanish, too. So we have 


\Pix{x),Pv{^)] = (A.6) 

For a free field and timelike x, = i{d^duA)/A. 

The above discussion can be extended to the case where is a collection of fields 

labeled by a (which could be spacetime indices), each conjugate to pa^- The obvious 

modification is: 

[r{x),Pb,iO)] = -iCl'^^%)d^A(x), (A.7) 

lPa/,(x),h^(0)] = -iC^lf^l(x)A(x). 

where S represents a product of Kronecker deltas. 
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